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Abstract 


An  equation  for  the  coalescence  of  colloidal  particles  under  the 
influence  of  velocity  gradients  has  been  applied  to  the  problem  of  coales¬ 
cence  of  spray  droplets  in  the  turbulent  fields  surrounding  two-fluid 

atomizers.  Two  derivations  have  been  made,  one  which  permits  the  estima¬ 
tion  of  the  probable  upper  limit  of  coalescence,  and  a  second  which  takes 
into  account  the  effects  of  dilution  of  the  spray  in  the  jet  which  reduces 
the  coalescence.  The  first  consideration  results  in  a  relationship  which 
may  be  evaluated  easily;  the  second  gives  rise  to  a  partial  differential 
equation  which  is  not  susceptible  to  analytical  solution.  A  solution  by 
iteration  of  the  corresponding  finite  difference  equation  is  discussed 
in  the  Appendix. 
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Introduction  *• 

^  q  ^  "•'  ,v> 

Th»  creation  and  use* of  fine0%prays  are  problem#  currently  receiving 

m  9  *•  * 

attention  in  several  fields  of  scientific  dndeavor.  Process  engineers  are 

concerned  with  systems  of  mass  transport  in  which  the  large  areas  involved 

>  * 

in  sprays  provide  a  means  for  increasing  overall  transfer  rates.  Engine 
designers,  for  similar  reasons,  pay  increasing  attention  to  fuel  sprays.  The 
stability  of  small  droplets  is  of  interest  in  military  problems  involving 
smoke  screens  and  the  dispersal  of  toxic  materials.  Agricultural  scientists 
are  similarly  concerned  with  dispersal  of  insecticides  and  chemical  sub¬ 
stances  which  control  or  regulate  the  growth  of  weeds.  In  all  of  these 
cases,  interest  in  the  use  of  sprays  leads  to  interest  in  methods  of  creating 
them. 

The  use  of  any  means  for  producing  a  fine  spray  demands  that  no  second¬ 
ary  effect  should  enlarge  the  spray  particles.  If  the  spray  is  produced  by 
a  two- fluid  nozzle,  there  is  little  to  be  gained  by  operating  the  nozzle  in 
such  a  way  as  to  form  extremely  fine  droplets  if  the  turbulence  generated  by 
the  nozzle  causes  the  droplets  to  coalesce.  Hence  the  problem  of  atomization 
by  a  two-fluid  nozzle  involves  both  the  design  of  the  nozzle  to  produce  spray 
of  the  desired  size,  and  the  prevention  of  coalescence  of  the  droplets  in  the 
turbulent  field  through  which  they  pass. 


Theory 

In  considering  the  effect  of  agitation  on  the  coalescence  of  colloidal 
systems  whose  particles  initially  are  spheres  of  uniform  diameter , 
Smoluchowski  (12)  derived  the  equation 

n  =  </VcO  (1) 

where  n  is  the  number  of  collisions  per  second  per  cubic  centimeter,  ^ 
is  the  concentration  of  particles  per  cubic  centimeter,  CO  is  the  velocity 
gradient,  centimeters  per  second  per  centimeter,  and  d  is  the  collision  diam¬ 
eter  of  the  particles  concerned,  in  centimeters.  The  collision  diameter  is 
defined  aB  the  distance  between  the  centers  of  two  spherical  particles  when 
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a  collision  occurs.  It  will  be  assumed  that  each  time  a  collision  occurs, 

coalescence  takes  place*  Since  each  collision  results  in  the  disappearance 

dj/ 


of  one  particle,  n  may  he  replaced  by  - 


d  t 


Furthermore,  as  a  result  of 


collision,  there  are  formed  particles  whose  collision  diameters  are  greater 
than  d,  and  such  particles  are  also  capable  of  removing  particles  from  the 
system  by  collision.  For  simplification,  it  will  be  assumed  that’ in  colli¬ 
sions  between  unlike  particles  the  collision  diameter  is  equal  to  the  diam¬ 
eter  of  the  larger  particle.  It  follows  that  the  rate  of  reduction  of 
by  collision  with  particles  of  all  sizes  is  given  by  the  equation 


Hi  ^3 


■h 


l 


-  4 
~  3 


CJV 


s  ■+* 


where  the  subscripts  refer  to  the  number  of  particles  of  the  original  size 
which  have  coalesced  to  form  each  particle  of  the  size  under  consideration. 
Now 

-jf-  [  4 +  »A  *  »Ai  ]=  f 

where  ^  is  the  loading,  or  the  volume  of  the  dispersed  phase  per  unit 
volume  of  mixture  of  both  phases;  therefore 


d  U1  zoo  Vj  nfr 


dt 


(2) 


The  nature  of  CO  depends  upon  the  hydrodynamics  of  the  system.  It  will  he 
assumed,  first,  that  the  velocities  of  the  dispersed  and  continuous  phases 
are  identical;  that  is,  no  slip  exists.  Tunitzki  (13)  defines  CJ  as  u/L 
where  u  and  L  are  the  turbulent  velocity  and  diameter,  respectively,  of  the 
smallest  eddy  capable  of  contributing  to  the  dissipation  of  energy.  Such  an 
eddy  is  characterized  by  a  lieynoids  number  of  unity  if  the  continuous  phase 
is  air: 


.  2 


hence 


-  CJ 
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where  ^  is  the  kinematic  viscosity  of  air.  Thus 
P 

M  _  en,J)  ni/juf-  (3) 

dt  ~  Ttu  r  i 

If  a  two-fluid,  nozzle  is  the  means  hy  which  the  spray  is  created,  the  system 
approximates  the  hydrodynamics  of  the  free  jet,  which  can  he  described 
adequately  hy  a  combination  of  theoretical  and  empirical  equations.  It  is 
first  necessary  to  assume  that  Smoluchowski ' s  equation  may  he  applied  without 
great  error  in  a  system  of  non-uniform  particles  to  the  particles  whose  diam¬ 
eters  are  within  an  arbitrarily  small  range  of  the  mean  diameter.  Equation 
(3)  will  therefore  be  used  without  subscript,  it  being  understood  that  the 
particles  under  consideration  axe  those  within  the  size  range  described 

above,  before  any  coalescence  occurs. 

Two  ways  of  treating  the  problem  now  appear.  In  one  we  may  attempt  to 
ascertain  an  approximate  upper  limit  to  the  coalescence  by  a  relatively 
simple  calculation.  As  an  alternative,  we  can  estimate  more  closely  the 
degree  of  coalescence  by  means  of  a  system  of  simultaneous  partial  differen¬ 
tial  equations  whose  solution  can  be  obtained  only  by  means  of  high-speed 
computing  devices.  Both  methods  will  be  investigated. 

Estimate  of  Upper  Limit  of  Coalescence 

The  upper  limit  of  coalescence  may  be  obtained  by  assuming  that  in 
Equation  (3)  the  term  Tf  is  constant  and  )J  is  affected  only  by  coales¬ 
cence,  and  by  applying  Equation  (3)  along  that  path  at  which  (u)  has  a 
maximal  value.  Such  conditions  exist  along  the  axis  of  the  jet,  if  the 
dilution  effect  of  the  jet  which  reduces  both  ^  and  id  is  ignored. 
Rather  than  consider  a  system  varying  with  time  as  Smoluchowski  did,  let  us 

consider  instead  one  varying  with  distance  along  the  axis  of  the  jet,  and  in 

.  dx 

a  steady  state  at  each  point.  This  may  be  done  by  expressing  u  as  ^  . 
Equation  (3)  then  becomes 

-  8VJL-  WlU) 


-  U 


(4) 
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Corrsin  (5)  has  shown  that,  along  the  axis  of  a  jet,  the  turbulent  velocity 

'll  rarely  exceeds  0.20  u.  Alexander,  Baron,  and  Comings  (2)  have  shown  that 
the  axial  velocity  u  of  a  free  jet  is  related  to  the  axial  displacement  x  by 
the  equation 


2 

u 


2  2 

n  c  x 


(5) 


where  Aq  is  the  area  of  the  jet  nozzle,  uq  is  the  velocity  at  the  nozzle, 
and  c  is  a  constant  whose  value  may  be  taken  as  0.075*  Substituting  Equa¬ 
tion  (5)  into  (4),  and  expressing  Aq  in  terms  of  the  diameter  of  a  round  jet, 

.tJL  .  (S> 

d  X  TT  n  2  c  X 

Separating  the  variables,  and  substituting  the  numerical  value  of  c, 


-  Hh =  0,68  T  uo B  ¥  w 

Since  the  turbulence  is  quite  low  in  the  potential  core  of  the  jet,  it  may  be 
assumed  that  no  coalescence  occurs  within  that  region;  hence  the  lower  limit 
of  x  should  be  4D,  where  experiments  show  the  potential  core  ends.  At  this 
point,  the  droplet  concentration  may  be  taken  as  u  %  the  original  value 
before  coalescence.  Then  Equation  (7)  may  be  integrated  to  give 


P 

In  — = — -  =  0.68  -  ''Iff  uqD  lnj 


Ji 


X 

4d  / 


Then,  recognizing  that  ucffi  is  Nfie  ,  the  Reynolds  number  for  flow  through 

H  6o 


the  nozzle,  the  above  equation  may  be  rewritten  as  follows: 


jj  ■  par 

\  / 


(s) 


Finally,  in  view  of  the  uncertainties  discussed  later,  which  are  inherent 
in  the  foregoing  derivation,  the  numerical  constant  0.68  should  be  replaced 


by  an  en^pirical  constant  I,  which  should  in  no  case  exceed  unity.  Thai, 


By  the  use  of  Equation  (8a)  the  snalleit  concentration  of  droplet*  likely 

to  exiit  along  the  axis  of  a  two-fluid  nossle  can  he  estimated,  provided  the 

loading  Off  ,  the  condition  of  airflow  at  the  noaile  as  characterised  by 

1L  ,  and  the  initial  concentration  of  droplets  of  the  aean  site,  lJ 0» 

*eo 

are  known. 

Correction  for  Dilution  in  the  Jet 

A  Closer  approxiaation  to  the  actual  aaount  of  coalescence  to  be 
expected  aust  take  into  account  the  reduction  in  caused  by  factors 

other  than  coalescence.  It  is  apparent  that  if  the  jet  entrains  air,  the 
concentration  of  droplets  aust  necessarily  fall,  and  the  loading  ^  aust 
also  be  reduced.  A  differential  ea^io“  »•*  UP  relating  the  rate  of 

accumulation  of  droplets  in  a  snail  volume  to  the  net  rate  of  transport  into 
the  volume  and  the  rate  of  generation  of  droplets  within  the  volumes 

4#-.  (9> 

<1* 

The  first  term  in  Equation  (9)  is  the  rate  of  accumulation  of  droplets,  the 
second  term  the  rate  of  transport  by  molecular  processes,  the  third  term  the 
rate  of  transport  by  convective  processes,  and  the  fourth  term  the  rate  of 
generation,  J0  is  an  effective  diffusivity  for  the  droplets  and  T  is  the 
velocity  vector.  At  steady  state,  the  rate  of  accumulation  is  sero,  and 
transport  by  molecular  processes  is  negligible  in  comparison  to  convective 
transports  thus  Iquation  (9)  becomes 

=0  (10) 

where  the  bar  indicates  the  time  average  of  the  quantity.  The  first  term  of 
Equation  (10)  represents  the  removal  of  droplets  from  the  element  of  volume 
by  convection,  and  the  second  the  creation  of  droplets  within  the  volume. 


Tor  the  cate  of  axial  symmetry,  tha  first  tera  above  nay  be  expanded  to  the 
fora 


V*Vt  -  ♦  i  (*  V  *>  (n) 

o  x  r  or 


Transport  perpendicular  to  the  axis,  expressed  in  the  the  third  tera  of 
Equation  (ll).  nay  now  be  treated  as  suggested  by  Baron  and  Alexander  (4), 

i  o  0  •  <f 


nr** 


C2x  d  17u 
2  dr 


(12) 


The  second  tern  of  Equation  (10)  is  the  rate  of  generation  of  droplets. 

Since  coalescence  removes  droplets,  the  rate  of  generation  is  negative.  This 
tera  represents  the  change  in  concentration  of  droplets  brought  about  by 
coalescence  alone,  and  is  given  by  Equation  (3),  which  will  be  rewritten  in 

the  fora 


Rj,  -  -  'f  (i)2  (13) 

Now  substitute  Equations  (ll),  (12),  and  (13)  into  (10),  there  results 
d  "J7u  Adi  d  17u\  _  8  i/p  qjf 

“FT  *  2  r  dr  |r  d  r  /  tr  »  ' 

The  fluctuating  velocity  u  nay  be  related  to  the  nean  velocity  u  at  any  point 
whose  co-ordinates  are  x,  r,  by  the  expression 

-2.  *  g(x,r) 
u 

vhtn  g(xpr)  is  fnpirio*!  function  of  x  wd  r0  Thoil 

(i>2  -  t2{®2  <15) 

Substituting  Equation  (15)  into  (l4) , 

|_  ©V  06) 

dx  2  r  d  r  I  dr  "  t* 


Assume  now  that  U  u  *  J^u 


(17) 


and  that  Ojf  u  *  ^u 
Equation  (l6)  way  then  he  rewritten 

VF-lHfV?)  <«> 

To  determine  at  any  value •  of  x  and  r,  Equation  (18)  muat  first  he 

solved  for  (  1/  u)  and  then  divided  hy  the  sui table  value  of  u. ~  This 
requires  a  knowledge  of  (  Hj/  u) ,  u  and  g  as  funotions  of  x  and  r. 

Equation  (18)  nil  obtained  hy  means  of  a  balance  over  the  number  of 
droplets  within  an  arbitrarily  small  sise  range.  A  similar  balance  dtn  be 
performed  with  respect  to  the  total  amount  of  the  dispersed  phase:  that  is, 

the  rate  of  accumulation  of  the  dispersed  phase  in  a  given  small  volume 
element  must  be  equal  to  the  net  rate  of  transport  of  that  phase  through  the 
surface  of  the  element  by  both  molecular  and  convective  transport  plus  the 
rate  of  creation  of  dispersed  phase  within  the  element.  The  mathematical 
statement  of  the  balance  is 

V  ‘Jj'vy  -  +  09) 

The  first  term  in  Equation  (19)  is  the  rate  of  accumulation  of  dispersed 
phase,  which  is  zero  in  a  steady-state  system;  the  second  term  is  the  rate 
of  transport  by  molecular  processes,,  which  should  be  negligible  in  comparison 
with  the  rate  of  transport  by  convective  processes,  the  third  term;  the  fourth 
term',  the  rate  of  creation  of  dispersed  phase,  is  zero.  Equation  (19) 
therefore  becomes 


V  •OjTI  «  0  (2°) 

If  Equation  (20)  is  expanded  for  the  case  of  axial  symmetry,  and  the  assump¬ 
tion  embodied  in  Equation  (12)  is  applied,  there  results 


10  - 


lunation  (21)  Is  linear  in  uj  Baron  and  Alexander  (4)  obtained  the 
following  equation  as  a  solution  at  all  values  of  r  and  x: 


In  Xquation  (22)  „  s  and  ^  refer  to  a  syste*  of  polar  co-ordinates  in  the 
plane  of  the  nozsle  fro*  which  the  jet  of  air  and  dispersed  phase  issue,  I)  is 
the  diameter  of  the  nozsle,  and  C  is  a  constant  whose  value  has  been  found 
empirically  to  be  0.087  (9)*  At  values  of  x  such  that  x/D  20 6  Equation 
(22)  simplifies  to  the  form 


By  similar  reasoning  applied  to  conservation  of  momentum,  there  results 


in  which  the  value  for  the  constant  c  has  been  found  to  be  0.0J5  (2).  At 
values  of  i  ^  20D,  Equation  (24)  simplifies  to  the  form 


By  the  same  assumption  as  was  involved  in  Xquation  (17),  u  may  be  obtained  by 

taking  the  square  root  of  the  result  of  Equations  (24)  or  (25). 

finally,  Cor r sin  (5)  has  measured  the  intensity  of  turbulence  at  various 

2 

points  in  the  field  of  a  free  jet.  from  his  data  values  off  g  have  been  com- 

o 

puted  as  a  function  of  x  and  r.  figures  1-3  show  g~  as  a  function  of  r,  at 
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e 


several  values  ©f  xj  figures  4  and  5  *how  t  aB  a  function  of  x  at  several 
values  of  r. 

Thus  jy  can  bo  determined  as  a  function  of  x  and  r  by  the  solution  of 
Equation  (18)  combined  with  Equations  (22)  and  (24)  and  the  data  of  figures 
1~5„  It  tes  not  been  possible  to  obtain  such  a  solution  analytically 0  and 

a  numerical  or  graphical  approach  is  indicated. 

The  methods  investigated  thus  far  are  integration  by  moans  of  the 
Differential  Analyser  (8),  by  punched-eard  computation  (11),  and  by  high¬ 
speed  digital  computer  (3).  The  two  last  mentioned  would  require  that  Equa¬ 
tion  (18)  be  set  up  as  a  finite-difference  equation.  This  has  been  done  by 
Hose  and  Sehilk  (ll)  for  this  case?  their  analysis  is  presented  in  the 
Appendix.  There  is  no  certainty  that  a  solution  can  be  obtained,  for  it  is 
by  no  means  obvious  that  the  difference  equations  converge  rapidly  enough  to 
give  a  solution  within  a  reasonable  number  of  iterative  computations  on 
punched-eard  computing  equipment.  Alt  (3)  bas  indicated  that  solution  br 
moans  of  a  high-speed  digital  computer  is  possible.  According  to  Hunt  (8), 
solution  by  means  of  the  Differential  Analyser  is  not  possible. 

Evaluation  of  Errors 

In  the  course  of  the  foregoing  derivations,  many  assumptions  have  been 
made,  some  of  which  obviously  are  not  completely  valid.  It  is  therefore 
necessary  to  investigate  the  effects  upon  the  resulting  solution  brought 
about  by  these  assumptions. 

The  equation  of  Smoluchowski  which  forms  the  basis  of  this  analysis  is 
subject  to  the  limitation  that  the  particles  concerned  be  of  such  a  site 
that  the  fluid  medium  behaves  as  a  continuum.  In  dealing  with  particles  in 
the  range  of  10  microns,  an  error  of  about  1  per  cent  is  made  by  ignoring 
the  particulate  nature  of  the  fluid?  this  error  Increases  markedly  as  the 
particle  sise  decreases,  and  decreases  as  the  particle  sise  increases.  The 
analysis  is  therefore  limited  to  particles  greater  than  10  microns  in 

diameter . 

The  assumption  of  spherical  particles  obviously  collapses  entirely  if 
the  particles  are  fibrous  or  needle-like.  If  the  particles  involved  are 
liquid  droplets,  the  error  introduced  by  the  assumption  of  sphericity  should 
not  be  great,  since  s^ll  liquid  droplets  are  often  spherical  or  nearly  so. 
There  is  also  question  about  the  assumption  that  coalescence  is  assured  each 
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tine  collision,  between  particles  occurs.  There  is  sons  doubt  that  this  is 
the  case  for  either  solid  particles  or  liquid  droplets.  On  the  other  hand , 
it  is  impossible  that  coalescence  occur  without  collision.  Thus  the  assump¬ 
tion  made  predicts  the  maximum  possible  rate  of  coalescence,  and  reduces  the 
probability  that  the  actual  coalescence  should  exceed  the  predicted  value. 

If  the  ratio  of  coalescence  to  collision  is  less  than  unity,  I  in  Iquation 
(8a)  should  be  less  than  0.68. 

The  assumption  that,  in  collision  between  unlike  particles,  the  colli¬ 
sion  diameter  is  equal  to  the  diameter  of  the  larger  particle,  is  obviously 
not  correct.  The  error  is  small  in  the  case  of  particles  not  widely  different 
in  size,  and  increases  with  the  difference  in  particle  size.  The  collision 
diameter  so  assumed  is  always  too  large,  and  the  effect  is,  again,  to  increase 
the  estimated  degree  of  coalescence.  However,  there  should  be  relatively  few 
particles  so  large  ttet  the  predicted  collision  diameter  is  seriously  in 
error.  Thus  the  number  of  collisions  involving  those  particles  should  be 
small.  In  any  case,  the  net  effect  is  to  bring  about  an  overly  large  calcu¬ 
lated  coalescence  rate,  which  is  safe  from  the  viewpoint  of  the  upper  limit 
calculation.  So  far  as  Iquation  (8a)  is  concerned,  this  would  make  the  true 

value  of  I  less  than  0.68. 

TuniUki's  relationships  regarding  the  term  CO  represent  the  most 
serious  possible  error.  Both  the  definition  of  CO  and  the  choice  of  a 
limiting  Reynolds  number  of  unity  are  without  substantial  experimental  confir¬ 
mation.  This  may  introduce  an  error  of  an  order  of  magnitude,  in  either 
direction. 

Finally,  the  application  of  the  Smoluehowski  equation  to  a  system  of 
particles  non-uniform  from  the  start  ignores  the  effects  of  coalescence  among 
particles  smaller  than  the  mean  diameter.  If  these  particles  coalesce  with 
each  other,  they  eventually  form  particles  of  the  mean  diameter,  and  bring 
about  a  decrease  in  the  total  rate  of  removal  of  particles  with  that  mean 
diameter.  If  they  coalesce  with  particles  larger  than  the  mean,  they  tend  to 
increase  the  diameters  of  those  particles,  and  hence  to  increase  the  rate  of 
removal  of  particles  of  the  mean  diameter.  It  is  likely  that  these  effects 
counterbalance  each  other  entirely,  or  that  the  net  effect  is  sc  small  as  to 
be  negligible. 
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In  summary,  the  assumption  regarding  the  relative  frequency  of  coales¬ 
cence  and  collision,  and  the  assumption  regarding  the  collision  diameter  of 
unlike  {articles  have  the  effect  of  increasing  the  estimated  degree  of 
coalescence  above  the  value  to  he  expected,  while  the  uncertainties  inherent 
in  the  value  of  CO  can  either  increase  or  decrease  the  estimated  value. 

It  appears  reasonable,  however,  to  suppose  that  all  these  possible  errors 
nay  be  encompassed  by  assuming  that  I  in  Equation  (8a)  nay  vary  froa  0.1  to 

1.0.  In  the  determination  of  the  upper  liait  of  coalescence,  the  assuap- 
tions  nake  the  systert  as  a  whole  behave  as  the  narrow  pencil  of  flow  down  the 
axis,  and  the  flowing  nature  of  the  jet  does  not  alter  the  loading  or  the 
droplet  concentration  within  the  pencil.  The  effect  of  this  choice  of  a 
aodel  assures  that  the  coalescence  predicted  will  exceed  that  which  occurs, 
and  that  the  upper  liait  predicted  will  be  greater  than  the  probable  value. 

The  correction  for  the  dilution  of  the  jet  eliminates  the  inaccuracy  of 
the  overly  siaplified  aodel  previously  described.  The  errors  inherent  in 
the  derivation  lie  in  the  extension  of  Corrsin's  data  to  yield  sufficiently 
detailed  intensity-of •turbulence  data,  in  the  averaging  procedures  by  means 
of  which  the  U  -  and  ^/-fluxes  transported  by  the  turbulent  components 
are  ignored,  leaving  only  the  fluxes  transported  by  aean  components  of  u, 

JlJ  ,  and  ,  and  in  the  rate  of  coalescence  term  as  expressed  in 

Equation  (13).  It  is  supposed  that  although  the  actual  amount  of  coalescence 
which  occurred  at  any  part  of  the  flow  field  may  be  subject  to  errors  as 
large  as  plus  or  minus  50  per  cent,  the  relative  amounts  of  coalescence  at 
various  parts  of  the  jet  should  be  reasonably  close  to  predicted  values. 

Thus  an  experimental  prograa  should  require  far  fewer  observations  than  would 
be  required  in  the  absence  of  a  guiding  relationship  such  as  that  here 
obtained. 

Discussion 

Ho  experimental  prograa  has  accompanied  this  work,  and  it  is  not  pos¬ 
sible  at  this  tine  to  obtain  a  quantitative  evaluation  of  the  relationships 
obtained.  There  are,  however,  experimental  fragments  which  are  in  qualita¬ 
tive  agreement  with  the  preceding  equations. 

It  lmts  been  observed  beyond  reasonable  doubt  that  coalescence  does  occur 
in  two-fluid  atomisation  systems.  Donnelly  (6)  has  shown,  by  means  of  high¬ 
speed  motion  pictures,  that  collision  and  consequent  coalescence  between 


single  droplets  in  &  spray  are  frequent  occurrences.  Hukiyama  and  Tanasawa 
(10),  in  their  equation  relating  mean  droplet  diameter  to  the  conditions  of 
operation  of  two-fluid  atomizing  nozzles,  show  that  the  particle  size 
increases  as  the  3/2  power  of  the  ratio  of  volumetric  liquid  flow  rate  to  the 
volumetric  gas  flow  rate.  This  is  consistent  with  the  predictions  here 
obtained;  one  might  therefore  postulate  that  the  increase  in  particle  size 
with  loading  is  a  phenomenon  related  to  coalescence. 

Teild  (j)  studied  the  charge  on  sprays  created  in  a  two-fluid  atomizer, 
using  cylindrical  probes  of  two  different  diameters,  connected  to  an  electron- 
ic  voltmeter ,  and  making  traverses  of  the  flow  field  at  several  different 
axial  positions.  He  integrated  the  total  charge  at  each  of  the  axial  posi¬ 
tions,  and  found  that  there  was  an  apparent  loss  of  charge  as  the  traverse 
moved  away  from  the  nozzle.  This  might  be  explained  a*  follows;  The  process 
of  atomization  wa®  accompanied  by  a  separation  of  charges,  thus  some  of  the 
droplets  retained  a  positive  charge,  some  a  negative  charge.  The  electronic 
voltmeter  was  sensitive  to  a  net  charge  of  one  sign  only.  If  all  the  drop¬ 
lets  in  existence  ©lose  to  the  nozzle  were  conserved,  then  there  should  have 
been  no  loss  of  charge,  and  successive  integrated  profiles  should  have  indi¬ 
cated  no  loss  of  charge.  If  coalescence  occurred,  it  is  probable  that  some 
droplets  bearing  opposite  charges  united,  and  that  neutralization  took  place. 
In  that  case,  the  resultant  droplets  must  have  had  a  net  charge  smaller  with 
respect  to  the  electronic  voltmeter  than  did  the  original  droplets  before 
coalescence.  Hence  the  loss  of  charge  is  consistent  with  the  occurrence  of 
coalescence.  Alexander  (l),  in  similar  experiments,  confirmed  the  observa¬ 
tions  noted  above. 

Xquation  (8a)  can  be  used  to  estimate  the  maximum  coalescence  likely  in 
a  given  application.  It  might  be  assumed,  for  example,  that  little  or  n© 
coalescence  would  occur  farther  than  thirty  nozzle  diameters  downstream  of 


the  nozzle.  Thus  JOD  might  be  substituted  for  x.  The  Reynolds  number  would 

be  determined  by  the  flow  conditions.  Then  a  particular  loading,  ' (JJ  » 

j  ' 

would  be  needed  in  order  to  prevent  the  ratio  — from  falling  below  the 

^ 0 

desired  value.  Or,  if  the  value  of  the  loading  were  predetermined,  it  would 

v  v  Jy 

be  possible  to  compute  the  smallest  value  of  the  ratio  .  likely  to  exist, 
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Tor  example,  suppose  that  it  was  necessary  to  sstiaats  tha  upper  Halt 
of  coalescence  for  a  ease  in  which  a  liquid ,  pr edoni nant ly  water,  was  atoa- 
ised  in  a  two-fluid  nossle  by  air  at  «0°f.  flowing  at  500  ft. /esc.  through  a 
l/2  inch  nossle.  The  gas  flow  rate  establishes  the  aean  droplet  else  orig¬ 
inally  foraod;  that  aean  sise  need  not  be  known  for  the  purposes  of  this 
calculation.  By  using  Iquatlon  (8a),  it  is  possible  to  coapute  the  fraotion 
of  the  droplets  of  the  original  aean  diaaeter  which  reaained  unaffected  by 
coalescence  at  a  distance  of  30  diameters  froa  the  nossle,  as  a  function  of 
yT  .  The  Reynolds  nuaber  for  the  flow  conditions  above  is  127,000; 
as suae  that  K  is  0.5»  A  weight  ratio,  water  to  air,  of  1,0  corresponds  to  a 
loading.  Off  ,  of  1.18  x  10"3„  Table  I  gives  as  a  function 

of  <y  •  0 


Table  I 

Upper  Liait  of  Coalescence;  jJ j jJ0 


as  a  function  of 

at  x/D  ■ 

r 

w 

1.18  x  104 

1.2  x  104^ 

1  x  104 

2.8  x  10'5 

1  x  10’5 

0.28 

1  x  104 

0.88 

The  solution  of  liquations  (18),  (22)  and  (24)  should  perait  one  to 
aake  a  close  eetiaate  of  the  actual  aaount  of  coalescence  likely  to  occur 
under  any  given  circuastances.  Because  of  the  complexities  of  the  solution, 
the  value  of  the  results  aust  be  weighed  against  the  expense  entailed.  Such 
an  evaluation  is  beyond  the  scope  of  this  report. 


Su—ary 

The  ecalescenoe  of  droplets  in  sprays  froa  two-fluid  nossles  has  been 
attacked  by  anally  with  the  coaleecence  of  colloidal  particlee  under  the 
influence  of  velocity  gradiente,  as  set  forth  by  Saoluchowskl .  Two  attacks 
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tove  been  employed:  the  first  considers  only  *  thin  paneil  of  spray  along 
the  axis  of  the  nozzle,  and  as  ranee  that  the  Jet  is  not  diluted  by  nixing 
with  air.  This  method  reeulte  in  an  upper  limit  to  the  amount  of  coalescence 
likely,  and  gives  riee  to  an  equation  which  may  be  readily  employed.  The 
second  method  takes  into  account  the  dilution  effect  of  the  Jet,  and  affords 
a  measure  of  the  number  of  particles  within  a  given  size  range  existing  at 
any  point  in  the  field  of  the  Jet.  The  latter  result  is  presented  in  the 
form  of  a  partial  differential  equation  whose  solution  by  analytical  means 
has  not  been  possible.  Solution  by  iteration  of  a  finite  difference 
equation  is  discussed  in  detail  in  the  Appendix. 

The  limits  of  accuracy  of  both  methods  are  discussed.  Experimental 
data  in  qualitative  agreement  with  this  treatment  are  described. 
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Table  of  nomenclature 


a 

Id 

Ao 

G 


c 


D 

XT 

J3' 


S 

h 

hr 

K 

L 


R, 


hr 

hr 

Area  of  a  round  nozzle 

Jet  spreading  coefficient  for  droplets,  whose  numerical  value 

is  0.08/ 

Jet  spreading  coefficient  for  momentum,  whose  nuraeripal  value 
is  0.0J5 

Diameter  of  a  round  nozzle 

Diffusion  coefficient  for  droplets  of  a  given  size 
Diffusion  coefficient  for  the  total  discontinuous  phase 
Diameter  of  particles  in  Smoluchowski ' s  equation  for  particle 
coale  sconce 

A  function  of  z  and  r  by  means  of  which  u  and  u  are  related 
Distance  between  two  points  on  a  grid 

Distance  between  two  points  on  a  grid  in  the  axial  direction 
Distance  between  two  points  on  a  grid  in  the  radial  direction 
A  constant 

A  length  characteristic  of  the  smallest  eddy  by  means  of  which 
energy  may  be  dissipated 

The  Reynolds  number  for  flow  through  the  nozzle 

The  rate  of  collision  per  unit  volume^  the  number  of  collisions 

per  second  per  cubic  centimeter 

The  rate  of  generation  of  droplets  of  a  given  size  per  cubic 
centimeter 

The  rate  of  generation  of  discontinuous  phase  per  cubic 


centimeter 


The  radial  and  axial  co-ordinates,  respectively ,  in  am  axially 
symmetrical  co-ordinate  system 

Normalized  radial  and  axial  co-ordinates  in  the  axially  symmet¬ 
rical  co-ordinate  system,  equal  to  r/S  and  x/D 
A  system  of  polar  co-ordinates  in  the  plane  of  the  nozzle 
Time,  seconds 

Axially  directed  velocity,  centimeters  per  second 
Fluctuating  component  of  axially  directed  velocity;  the  root- 
mean-square  value  of  the  velocity  fluctuation 
The  velocity  vector 

The  viscosity  of  air,  grams  per  centimeter  per  second 
The  concentration  of  droplets  of  a  given  size,  in  droplets 
per  cubic  centimeter 

The  density  of  air,  grams  per  cubic  centimeter 
The  loading,  volume  of  the  dispersed  phase  per  unit  volume  of 
mixture  of  both  phases 
A  velocity  gradient,  in  sec.1 

mv  .  *  d(  )  ,  Id  (r  )  .  1  d  (  ) 

The  vector  operator,  5*  +  "t~Kt  +  r  *d6 
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Appendix 


BEPOHT  ON  FHBLIMINAHY  INVESTIGATION  OF  SOLUTION  Of 
ILLINOIS  JET  COALESCENCE  EQUATION 

Arthur  Hose  and  Joan  A.  Schilk 
The  Pennsylvania  State  College 
State  College,  Fa. 

The  method  of  iteration  has  been  carefully  studied  as  a  possible  approach 
to  the  solution  of  the  jet  coalescence  differential  equation 


c)j>U 


c%  d  f l -i  dm 


2.  r 


dr 


d  m]  _ 
dr 


&£_ 
If  *4 


(A-l) 


This  iteration  method  involves  first  the  estimation  of  a  series  of  first 
trial  values  of  the  dependent  variabl e ,  and  then  improvement  of  these  by 
means  of  successive  iterations  with  the  difference  equation  corresponding  to 
the  differential  equation. 

If  j“  =  (IH),  and  f(x»r)  *  (  r\j/'  u)g^»  then  Equation  (A-l), 
expressed  in  terms  of  difference  quotients,  takes  the  following  forms:  Using 

the  average  first  difference  quotients,  and  hr  =  hr 
(see  nomenclature) 

when  r  J  0 


(%r) 


d%h 

C  -t- b($b'htr)- rfy r)j 

l  ■+  #5?  fit  r ) 

™  (M) 


when  r  =  0 


2.CVAi*  r ->-h) 


+  g  r  -k, 


'J]+  ')  -  §(#  -h  r  3 


4 C~%  +  4{7r) 


(A-3) 


I  h>r)  = 
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when  r  =  0 


In  order  that  the  difference  equation  can  he  used  for  solution  to 
the  problem,  a  grid  with  initial  trial  values  of  the  dependent  variable, 

^  ,  must  be  constructed.  The  iterative  method  of  solution  merely 

Involves  repeated  use  of  the  difference  equations  for  one  point  after  another 
in  the  grid,  and  repetition  of  this  procedure,  until  successive  values  of 
j°(/Xtr7)  for  a  given  point  are  as  nearly  identical  as  is  desired. 

For  the  present  problem,  a  grid  for  integral  values  of  'p  =  0  to  40  and 
r  =  0  to  5  was  necessary  because  of  the  boundary  conditions  (see  Fig.  6). 
This  grid  was  first  prepared  for  hjt.  -  hp  =1.0,  obtaining  the 
first  trial  values  of  by  linear  interpolation  from  the 

boundary  values  (see  Table  II).  Values  of  f(x,r)  were  also  calculated  for 
the  points  of  this  grid  (see  Table  III).  Using  Equations  (A-2)  and  (A-3) 
and  h  »  1.0,  a  series  of  pilot  calculations  were  made  at  various  points 
of  the  grid.  This  involved  two  iterations  and  use  of  Equations  (A-2)  or 
(A=>3 ) ,  a  total  of  115  times  to  obtain  115  aew  values  of  <k,  .It  was  found 
that  some  of  the  resulting  values  of  were  greater  than 

1.0,  or  were  negative.  In  an  attempt  to  avoid  this  lack  of  convergence,  the 
size  of  h  was  reduced.  Additional  first  trial  values  of  £ ('ft r) 
were  obtained  for  fractional  values  of  F  and  %  »  by  linear  interpola¬ 

tion  from  previous  values,  and  more  pilot  calculations  were  then  made. 
Reduction  in  h  greatly  increases  the  number  of  points  on  the  grid,  and 
the  amount  of  computation  required. 

The  following  values  of  b)  ,  'P  .and  P  were  used  for  the  succes¬ 
sive  series  cf  trial  calculations,  with  the  objective  of  finding  an  approach 
that  would  lead  to  convergence  (see  Table  IV). 

These  pilot  calculations  showed  that  the  largest  value  of  h  which 
could  be  used  in  Equations  (A-2)  and  (A=>3)  was  ^  =0.1.  This  would 
involve  calculation  of  for  approximately  12,000  points  for 


Using  forward  first  difference  quotients  and 
when  r  4  0 
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Using  backward  first  difference  quotients  and  /7^ 
when  H  4  ® 
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Using  average  first  difference  quotients, 

hr=  b 

when  P  4  0 
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each  iteration.  This  ie  a  prohibitive  volume  of  computation. 

Use  of  Equations  (A— 4)  and  (A-5)  and  (A- 6)  and  (A-7)  resulted  in 
reasonable  values  of  £  for  b)  =  1.0  in  most  portions  of  the 

grid,  but  resulted  in  nonconvergence  near  certain  of  the  boundaries. 

Trial  calculations  were  next  performed  with  =1.0  and 

hr*  -  0.1.  It  was  found  that  this  method  resulted  in  convergent 
values  of  throughout  the  grid,  as  far  as  the  process  was 

carried. 

A  tentative  IBM  program  plan  has  therefore  been  worked  out  on  this 
basis,  and  cost  estimates  prepared  accordingly. 

As  far  as  can  be  determined,  liquation  (A-l)  is  a  type  of  partial 
differential  equation  whose  numerical  solution  has  not  been  attempted 
previously.  There  is,  therefore,  no  guarantee  that  the  proposed  iterative 
computation  will  give  convergence  and  a  solution  within  a  reasonable  number 
of  iterations.  Apparently,  the  solution  of  this  equation  is  itself  a 
research  problem  of  considerable  importance  in  applied  mathematics. 
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TABLE  II 

ABBREVIATED  TABLE  0?  FIRST  TRIAL  VALUES  OF  £(*  r) 


%/r 

0 

1 

2 

3 

4 

5 

i 

•9750 

0 

0 

0 

0 

0 

5 

.8750 

.2500 

0 

0 

0 

0 

10 

15 

■m 

.2813 

.3125 

0 

0 

0 

0 

0 

0 

0 

0 

20 

.5000 

•3077 

.1154 

0 

0 

0 

25 

.3750 

.2613 

.1477 

.0341 

0 

0 

30 

.2500 

.1859 

.1218 

•0577 

0 

0 

35 

.1250 

.0979 

.0707 

.0435 

.0163 

0 

39 

.0250 

.0201 

.0152 

.0103 

.0054 

.0005 

Notes  The  complete  table  (12J  values  of  £  )  has  been 

prepared  for  all  intermediate  integral  values  of  X 
and  r  ,  and  also  for  certain  additional  intermediate 
fractional  values. 


TABLE  III 

ABBREVIATED  TABLE  OF  FIRST  TRIAL  VALUES  OF  /(ft/*)  x  10^ 


=  1  x  10 


%jr_ 

0 

1 

2 

3 

4 

.,.,5 

1 

0 

0 

0 

0 

0 

0 

5 

1.515 

0.053 

0 

0 

0 

0 

10 

1.427 

0.234 

0 

0 

0 

0 

15 

0.832 

0.541 

0 

0 

0 

0 

20 

0.470 

0.388 

0.153 

0 

0 

0 

25 

0.302 

0.271 

0.155 

0.054 

0 

0 

30 

0.210 

o.iq6 

0.137 

0.068 

0 

0 

35 

0.154 

0.149 

0.116 

0.071 

0.033 

0 

39 

0.124 

0.122 

0.101 

0.068 

0.038 

0.017 

Notes  The  complete  table  (12/  values)  has  been  prepared  for  all 
intermediate  integral  values  of  %  and  r  ,  and  also 
for  certain  additional  intermediate  fractional  values. 
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Figure  4.  Turbulence  Function  vs.  Normalized  Axial  Co-ordinate 
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